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ABSTRACT 



Context. Weak entangled magnetic fields with mixed polarity occupy the main part of the quiet Sun. The Zeeman effect diagnostics 
fails to measure such fields because of cancellation in circular polarization. However, the Hanle effect diagnostics, accessible through 
the second solar spectrum, provides us with a very sensitive tool for studying the distribution of weak magnetic fields on the Sun. 
Aims. Molecular lines are very strong and even dominate in some regions of the second solar spectrum. The CN B 2 X - X 2 X system is 
one of the richest and most promising systems for molecular diagnostics and well suited for the application of the differential Hanle 
effect method. The aim is to interpret observations of the CN S 2 Z - X 2 E system using the Hanle effect and to obtain an estimation of 
the magnetic field strength. 

Methods. We assume that the CN molecular layer is situated above the region where the continuum radiation is formed and employ 
the single-scattering approximation. Together with the Hanle effect theory this provides us with a model that can diagnose turbulent 
magnetic fields. 

Results. We have succeeded in fitting modeled CN lines in several regions of the second solar spectrum to observations and obtained 
a magnetic field strength in the range from 10-30 G in the upper solar photosphere depending on the considered lines. 

Key words. Line: formation - Sun: magnetic fields - Molecular processes - Polarization - Radiative transfer - Scattering 



^ 1. Introduction 

■ The understanding of the solar surface magnetic field plays a 
| key role in modern solar physics since the magnetic field is 
, connected to and even controls most phenomena of solar ac- 

f^) ■ tivity. The solar magnetic field is structured in a fractal-like 
t***'- \ pattern with a significant degree o f self-similarity over many 
O ■ scales dStenflo & Holzreuteil l2003l) . while the field strength 
^ covers many orders of magnitude reaching up to a few kG. 
• i-h , Observationally, we can resolve only a small fraction of the 

■ magnetic structures, which can be accurately measured via the 
5h ' Zeeman effect diagnostics. However, about 99% of the photo- 

. 5^ i spheric volume is filled by unresolved, entangled magnetic fields 
with mixed polarity, not accessible with the Zeeman effect due 
to cancellation of the polarimetric signal, that might contribute a 
substantial amount of ma gnetic energy to the solar photosphere 
dTruiillo Bueno et alj|2004 . 

Th ese hidden ma gnetic fields can be accessed with the Hanle 
effect dStenfloiri982l) . which modifies the properties of coher- 
ent scattering processes. The Hanle effect can thus only be em- 
ployed through the spectrum formed by the c oherent scatte ring, 
which is called the "second solar spectrum" dlvanovlll99ll) due 
to its richness and significance. The Hanle effect usually leads 
to a reduction and rotation of line polarization as compared 
to the non-magnetic case. The net-influence of rotation cancels 
out in a turbulent magnetic field, but the depolarization remains 
visible for most geometries even in a turbulent magnetic field. 
Moreover, the Hanle effect is extremely sensitive to the weak 
field. Therefore, the second solar spectrum provides us with a 
very sensitive tool for studying the distribution of weak magnetic 
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fields on the Sun and gives the possibility to measure spatially 
unresolved mixed polarity magnetic fields. 

Molecular lines are particularly valuable for magnetic field 
diagnostics thanks to their broad range of magnetic sensitiv- 
ities within narrow spectral regio ns. This allows us to em- 
ploy the differential Hanle effect dStenflo et all 1 1998b . which 
greatly reduces the model dependence of deduced magnetic 
field strengths, and has lead to the first unambiguous de- 
tectio n of molecular Hanle eff ect in the second solar spec- 
trum (Fluri & Berdyugina 2004). Recent studies of the molec- 
ular Hanle effect have yielded field strengths of the order 
of 10 G for the quiet solar photosphere, assuming a single- 
valued, isotropic, turbulent magnetic field (F aurobert & Arnaud 
120031: iBerdvugina & Flurill2004t iTruiillo Bueno et alj|2004l) . al- 
though larger field strengths are possible depending on the as- 
sumed collision rate, in particular when interpreting MgH line s 
dAsensio Ramos & Truiillo Buenoll2005l:lBommier et al.ll2006T) . 

Diagnostics with molecular Hanle effect has concentrated so 
far only on two molecules, C2 and MgH. It is very important to 
extend the technique to other molecules because in general dif- 
ferent molecules form at different heights due to their large tem- 
perature sensitivity and can thus sample other layers of the solar 
atmosphere. Combined with atomic Hanle effect, wh ich is well 
established since the work by Faurobert -Scholll (fT992), only such 
a multitude of lines will allow us to constrain the 3 -dimensional 
structure of the unresolved solar magnetic field. 

In the present paper we apply the Hanle effect to a new 
molecule, namely the CN violet system, for which we have re- 
cently developed a complete theory of the Hanle effect taking 
into account the Pasc hen-Back effect and interference between 
doublet components (Shapiro et al. 2007). After describing our 
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radiative transfer method we fit our model to observations in 
three spectral regions of the second solar spectrum, which en- 
ables us to determine the magnetic field strength in the quiet 
solar photosphere. 



2. Radiative transfer model 

In this section we introduce our method for solving the radiative 
tra nsfer problem. It is based o n the modeling strategy developed 
by iFaurobert & Arnaudl (120031) . who have calculated scattering 
polarization at the line centers of C2 and MgH transitions. Here, 
we expand their method to obtain wavelength dependent spectra 
not restricted to line centers. 

The radiative transfer model is optimized for relatively 
weak lines whose scattering polarization is calculated within 
the framework of a plane-parallel atmosphere using the single- 
scattering approximation. The CN lines are assumed to form in 
a homogeneous atmospheric layer situated above the continuum 
formation region. We fully account for line blending and for de- 
polarization of the continuum polarization by line absorption 
and scattering. In this respect, t he current model is more com- 
plex than the one employed by Berdvugina& Fluri (2004]) for 
very weak C2 lines. The CN bands that are modeled in this pa- 
per include stronger lines and require a consistent treatment of 
continuum depolarization, although we increase at the same time 
the model dependence of the magnetic field diagnostics, because 
we rely on a proper computation of the continuum polarization. 
Nonetheless, the model remains simple and ideally suited for 
differential Hanle effect diagnostics. 

First we introduce the necessary notation and the source 
function. Then we discuss in detail the two-layer approximation, 
the free model parameters, and the fitting procedure for compar- 
ing the model with observations. 

2.1. Source function 

The total, monochromatic optical depth is given by 



dr(v) = - 



Ur + 



J>y'(v) 



dz, 



(1) 



where S = S(t, v, jj) is the total source function 



k c B th + I ky(v)Bto + cr c S c + £ cr^ (v)S* 

/ = 1 

n _ , 



(4) 



with B t h the Planck function, S c =S c {t, v,fi) the continuum source 
function, and S' e =S' ( (t, v,/j) the scattering part of the line source 
function of the z'-th line. 

In our further discussion we employ the single-scattering ap- 
proximation, since the CN lines for which we will apply this the- 
ory are optically thin. Then, the scattering part of the line source 
function takes the form 

1 

Sj = \ J Ffa//, % n)/i„cOu') exp(- Tmc(v ?~ T ) dyt, (5) 


where /i n c=J(Ti n c, v, fi) is the Stokes vector at the lower boundary 
with an optical depth Ti nc . 

Within the weak-field approximation and under the assump- 
tion of an isotropic turbulent magnetic field the phase matrix P' ( 
is parametrized as 



r t = £n + -w' 2 w{ l p (2) . 



(6) 



The matrix En represents isotropic, unpolarized scattering, 
while P (2) des cribes linearly polarized coherent scattering (e.g. 
j Stenflol 1 19941). T he line polarizability W' 2 , calculated as in 
Berdvugi na et al.l (|2002), corresponds to the branching ratio of 
polarized and unpolarized scattering. The physics of the Hanle 
effect is contained in the Hanle depolarization factor 



Wn = 1 - 0.4 



with 



y H = 0.88 



8lB 



r R + D' 2 ) ' 



(7) 



(8) 



where a c — k c + cr c is the continuum opacity (with k c and <x c the 
continuum absorption and scattering coefficients), a' ( = k' t + ai 
is the integrated line opacity in the z'-th line (with k' t and <x' f the 
line absorption and scattering coefficients), and <p'(v) is the nor- 
malized Voigt profile function. The opacities in molecular lines 
depend on the main line parameters and on the temperature of 
the CN layer T eS 



4 K g'ef'e exp(-£y£r eff ), 



(2) 



where g'. and E' ( are the statistical weight and energy, respec- 
tively, of the lower level, and fl is the absorption oscillator 
strength, wh ich can be calculated with the pure Hund's case (b) 
expressions dSchadedll97llWhiting & Nichollslll974 because 
the orbital angular momenta for both the upper and lower states 
are equal to zero in the CN violet system. 

Now we can write the radiative transfer equation in the usual 
form 



=I~S, 

OT 



(3) 



The variables Fr and D (2) are the radiative damping and depo- 
larizing collision rates, respectively, and B is the strength of the 
single-valued, isotropic, microturbulent magnetic field. For the 
analysis in this paper we neglect depolarizing collisions. In the 
CN violet system we are dealing with 2 2 - 2 E transitions, so that 
the effective Lande factors gL of the upper state of the two dou- 
blet components can be expressed through the quantum number 
of the total angular momentum J or, alternatively, of the angular 
momentum N, i.e. the total angular momentum without spin, 



gu = \I{N+\) = \U, 

fta = -!/(# +5) = -!/(/+ 1), 



(9) 



where the indices 1 and 2 refer to the first and second doublet 
components with J - N- 1 /2 and with J - N+ 1 /2, respectively. 

With our choice of the phase matrix we implicitly neglect in- 
terferences between different fine structure components and per- 
turbations by the Paschen-Back effect. The latter becomes rel- 
evant in th e CN violet system only for field strengths exceed- 
ing 100 G dShapiro et al.ll2007l) . while the fields detected in the 
present study are weaker than 30 G. 
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Let us rewrite Eq. © in the form 

n 

0"c + 2 CT>'(V) ( „ 



a c + 2 a' ( <p'(v) 



\ (=i 



+ — flth, 

« in? 



where the coefficients Xi an d Xc correspond to the conditional 
probabilities that scattering occurs in the j-th line or in the con- 
tinuum, respectively. In general they depend on frequency and 
are defined as 



Xi = 



a-' c ip'(v) 



Xc = 



(11) 



cr c + £ CT J e (pJ(v) 



The emergent Stokes vector is given by the formal solution of 
the radiative transfer equation 



TW(V) 



/ em (v,ju)= f S(r,v,/i)e^— +/ inc (v, J u)e- T ^. (12) 
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Fig. 2. The dependence of the Hanle depolarization factor Wh, 
given in Eq. @, on the magnetic field strength and the N quan- 
tum number for the CN violet (0,0) band. The three curves cor- 
respond to the depolarization factors Wh equal to 0.9 (solid), 0.5 
(dashed), and 0.3 (dotted), the last of which is already close to 
the asymptotic saturation limit of 0.2. 



2.2. Two-layer approximation 

We assume that the CN molecular layer is situated above the re- 
gion where the continuum radiation is formed. The CN layer 
is isothermal with a homogeneous CN density, and its lower 
boundary is irradiated by the solar continuum radiation. While 
passing through the molecular layer this initial incident radiation 
can be scattered (but only once, as we adopt the single scattering 
approximation) or absorbed by CN molecules. 

The "continuum" layer corresponds to the lower photosphere 
of the average quiet Sun, defines the lower boundary condi- 
tion of the CN layer, and is thus crucial for the formation of 
the CN lines. In particular, the continuum layer determines the 
anisotropy of the continuum radiation, without which scattering 
polarization would cancel out, and it sets the degree of contin- 
uum polarization, which is subsequently depolarized by the line 
opacity. 

In our modeling strategy we take the continuum intensity 
and polarization from previous work, rather than recalculating 
these values. The continuum intensity is obtained from the an- 
alytical expressions fo r limb-darken ing derived from observa- 
tions of the quiet Sun (lNeckellll996h . We assume the radiation 
field to be azimuthally symmetric. The value s of the conti nuum 
polarization correspond to the calculations bv lStenflol (120051) for 
jj = 0.1. The continuum polarization for arbitrary //-values are 
computed by scaling the value at u — 0.1 with t he analytical ex- 
pressions introduced by iFluri &Stenflol dl999). which describe 
the yu-dependence of the continuum polarization and are based 
on theoretical modeling of the quiet solar photosphere. 

The CN layer is described by the source function and the 
radiative transfer equation given in Sect. 12.11 In Eqs. (0 and 
(fT2l we insert the Stokes vector /i nc , which is incident into the 
CN layer at the lower boundary and given by the frequency and 
angle dependent continuum intensity and polarization formed in 
the continuum layer. The same equations require also the total 
optical thickness Ti nc = 7~cn of the CN layer. Assuming a c 
0, since the continuum is formed below the CN layer, the total 



optical thickness of the CN layer Ten can be calculated with the 
expression 



TcN (V) = Tnoim ■ ^ a' ( (f>'(v) 
1=1 



(13) 



where the frequency independent coefficient 7~ norm is a free pa- 
rameter of the model. It depends on the CN density and the geo- 
metrical thickness of the CN layer. 

Within the two-layer approximation the total source func- 
tion, Eq. ( TT7JI ). is further simplified since the opacity in the CN 
layer is only due to CN molecules, i.e. k c = cr c = 0. We as- 
sume that the ratio of the scattering and absorption probabilities 
c sc =a i/ a t is me same for all molecular lines. This is basically 
equivalent to the assumption that the inelastic collision rate is 
independent of the involved transitions because the life time of 
the exited CN state B 2 1. is almost independent of its quantum 
numbers (with the exception of some states with very small J 
numbers). Furthermore we neglect in our model the elastic col- 
lisions. Equation dTOb is then reduced to the form 



V i=l 



+ (l-c sc )B, 



th • 



(14) 



Here the first term with the summation represents the scattered 
diffuse part of the radiation field, restricted within the CN layer 
to line contributions, while the second term corresponds to the 
thermal part. 

The emergent Stokes vector is calculated employing the for- 
mal solution of the radiative transfer equation. To improve the 
accuracy of our results we divide the CN layer into k sublay- 
ers with equal optical thickness. The source function is assumed 
to be constant within each sublayer and evaluated at the center 
of each layer, giving us k different values Si. . .S*. Introducing 
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30 40 
Angular momentum N 

Fig. 1. Differential Hanle effect in synthetic Stokes profiles of R-branch lines in the CN violet (0,0) band. For selected angular 
momenta N of the lower state we show the complete doublet, which causes the double peak structure for large N. Three values for 
the turbulent magnetic field strength were chosen, giving identical Stokes / profiles but large modifications in the Stokes Q profiles: 
0G (upper curve), 10 G (intermediate curve), and 100 G (lower curve). In this and all other figures we assume an elastic collision 
rate equal to zero. 



the above assumptions into Eq. ( [TZl l. we find for the emergent 
radiation field at a given frequency 



IQi) = /incOu)e 



■Ten If 



) -(15) 



7=1 



The first term describes the radiation transmitted through the 
CN layer while the summation in the second term adds up the 
contributions from the different CN sublayers to the emergent 
spectrum. Note that the depth dependence of the source func- 
tion within the CN layer is introduced by the exponential factor 
in Eq. (0, which describes absorption of the incident radiation 
and thus the height variation of the radiation field even within 
the isothermal slab. The scattering coefficient, on the other hand, 
depends only on the probability of the collisional de-excitation 
of the upper CN states, which under our assumptions remains 
constant for the whole molecular layer. 

The subdivision of the molecular layer, expressed by 
Eq. (|T3T>, is a further sophistication of the model introduced by 
Faurobert & Arnaud (2003), who have assumed one single value 
of the source function within the molecular layer. The sublayers 
take care of the depth dependence of the source function within 
the CN layer, which results from the optical depth r in the ex- 
ponent of Eq. (0. Note that this does not contradict to our as- 
sumption of an isothermal CN layer with constant density. Our 
test calculations have shown that the emergent Stokes vector be- 
comes independent of the number of sublayers for k greater than 
five. 



2.3. Free parameters 

The model described above contains several free parameters, 
which have to be determined by fitting the theoretical spectrum 
to observations. These are the temperature of the isothermal CN 
layer Teg, the scattering probability c sc , the optical depth nor- 
malization coefficient 7^ rm> an d the magnetic field strength B. 

The temperature affects the emergent spectrum in two 
ways. First of all, it enters via thermal term B t h of the source 
function (cf. Eq. (TBI). This strongly influences the intensity, but 
not the other Stokes parameters, because of the strong temper- 
ature dependence of the Planck function and because the con- 
tribution from the thermal part is relatively large in the case of 
CN lines. Secondly, the line opacities depend on temperature via 
the Boltzmann factor (cf. Eq.|2]l, which affects the whole Stokes 
vector, although in a much weaker way than the thermal term. 

The scattering probability c sc defines the balance between 
the thermal and nonthermal parts of the total source function 
and influences the whole Stokes vector. The contribution to the 
Stokes vector from the CN layer linearly depends on c sc . 

The optical depth normalization coefficient 7~ n0lm describes 
the relative contributions to the Stokes vector from the CN layer 
and the continuum layer. It effectively determines the optical 
thickness of the CN layer, thus scaling both the molecular line 
depths in the intensity spectrum and the degree of line polariza- 
tion. 

Our main goal is the determination of the magnetic field 
strength B assuming a turbulent, isotropic field distribution 
within the CN layer. The magnetic field enters via Eq. © into the 
line source function. Since we consider weak fields and neglect 
Zeeman broadening, the magnetic field strength has virtually no 
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effect on Stokes / and only causes differential Hanle effect in the 
second solar spectrum (see Sec. 13. Q , 

Usually, the influence of B on the Q/I c spectrum can be eas- 
ily distinguished from the effects of c sc and T noTm . The turbu- 
lent magnetic field causes a depolarization of the CN lines. The 
same can in principle also be obtained with a smaller scattering 
probability c sc . Even an increased T noim parameter can mimic 
the same behavior because it leads to stronger continuum depo- 
larization. The crucial point, however, is that c sc and T noim affect 
all molecular lines with comparable strength, while the magnetic 
field leads to differential depolarization, because the Hanle effect 
efficiency strongly depends on the J quantum number (Fig. [T]). 
Thus, the relative strength of CN lines in the second solar spec- 
trum mainly depends on the magnetic field, but not so much on 
the other free parameters. Therefore, it is possible to diagnose 
the magnetic field strength when fitting several lines with differ- 
ent Hanle effect sensitivities. 



2.4. Fitting procedure 

We simultaneously fit the modeled polarization spectrum to ob- 
served Q/I c data and the modeled intensity spectrum to /// c 
observations. During the fitting procedure we vary the mag- 
netic field strength B, the scattering probability c sc , the effective 
temperature T e s and the optical depth normalization coefficient 
T noim . The most suitable values for these free parameters are 
determined with a least-square method. More precisely, we min- 
imize the quantity 



X 



2 _ i v/(e/«ob S -(e//c)ti 



= -Y 



°7//c 



(16) 



where n is the number of data points and o-q/ Ic and 07^ are the 
absolute errors of the observed polarization and intensity, which 
we set to 0.008% and 0.005%, respectively. The coefficient (, 
with possible values in the range from to 1, allows us to give 
larger weight either to the intensity or polarization fit. As the sig- 
nal in polarization is smaller and as only polarization is affected 
by the magnetic field we overweight the fit to the polarization fit 
by choosing £ = 7/8 (but the result of the fitting procedure is 
almost independent on the assumed ( value). 

The x 1 value in Eq. ( fTrjj l is a measure of the discrepancy be- 
tween observations and the theoretical model, and describes the 
accuracy of the fit. Wavelength points that are heavily affected by 
atomic blends are excluded from the summation and the fitting 
procedure since we neglect atomic line opacity in our model. 



3. Diagnostics with the CN violet system 

In this section we apply the previously developed theory to the 
interpretation of observations of the CN B 2 E - X 2 !. system, the 
so-called violet system, in several spectral regions, with the goal 
to determine the turbulent magnetic field strength in the quiet 
Sun. From the discussion in Sect. 12.31 it is clear that suitable 
wavelength regions should contain molecular lines with differ- 
ent Hanle effect sensitivities. Therefore, we first investigate the 
differential Hanle effect in the CN violet system. This will guide 
us to choose three ideal spectral windows in different molecular 
bands, which are then analyzed and discussed. 
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Fig. 4. Fits to the observations (thin solid lines) of Stokes I/I c 
and Q/I c in the first wavelength region (3865.9-3866.8 A). The 
assumed magnetic field strengths in the calculated spectra are 
G (no depolarization, thick dashed line), 27 G (best fit, thick 
solid line), and 45 G (too much depolarization, dotted line). The 
labels at the bottom of the lower panel indicate the position of 
the lines in the CN doublets with the corresponding N quantum 
number of the lower states. The labels at the top of the lower 
panel indicate the wavelength points which were chosen for the 
fitting procedure. 



3. 1 . Differential Hanle effect 

Figure[TJillustrates the different Hanle effect sensitivities of syn- 
thetic line profiles for the example of the CN violet (0,0) band 
(other CN violet bands give the same qualitative results). The 
relative degree of line polarization in the non-magnetic case 
is defined by three factors: the oscillator strength (which in- 
creases with the J number), the scattering polarizability W2 
(which decreases with J and asymptotically approaches a value 
of 0.1), and the Boltzmann factor (which also decreases with 
J). However, the final polarization obtained within the lines not 
only depends on the intrinsic line polarization. Rather it results 
from a combination of different contributions, namely the con- 
tinuum polarization, continuum depolarization due to the line 
opacity, and line polarization, where only the latter is subject to 
the Hanle effect. It is clearly visible in Fig. [T]tha t stronger lines 
lead t o additional continuum depolarization (cf. iFluri & St enflo 
120031) . Furthermore, we can identify the expected diverse behav- 
ior of the lines in a magnetic field. While the intensity profiles 
remain unaffected, the line polarization is depolarized with in- 
creasing field strength, reflecting the J dependence of Wu- A 
magnetic field of 10 G already leads to a nearly saturated Hanle 
effect in the lines with N = 1 and 5, while the polarization in 
the N = 70 and 80 lines remains almost unchanged compared to 
the non-magnetic case. The latter lines require a magnetic field 
strength in the order of 100 G for a considerable depolarization. 

The magnetic field necessary for sizable depolarization by 
the Hanle effect is shown in Fig.|2]for different N quantum num- 
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Fig. 3-x 2 - contours of the fitting procedure in the first wavelength region (3865.9-3866.8 A). The panels show 2-dimensional cuts 
in the parameter space through the minimum^- 2 . = 1.13. Numbers on the contours indicate the values of^ 2 as defined in Eq. (fT~6b . 



bers, again for the CN violet (0,0) band. The different curves 
correspond to three values of Wn and identify the magnetic field 
range in which a specific line is sensitive to the Hanle effect. 
For example, lines with N = 40 react to the Hanle effect mainly 
in the field range from about 10 G to 60 G. For stronger fields 
these lines are maximally depolarized. On the other hand, their 
Hanle effect depolarization is negligible for weaker fields, which 
makes them ideal reference lines below 10 G when compared to 
lines with much smaller N. 

The ideal choice of spectral windows thus depends on the 
actual magnetic field strength. Based on previous results we 
selected narrow wavelength regions containing lines with N 
greater than 40 and others with N around 10. An additional se- 
lection criterion is a small number of atomic or molecular blends 
since they would significantly alter the intensity and polarization 
spectra. 

3.2. Observational data 

We compare our model calculations with the data f r om th e 
second solar spectrum atlas compiled by iGandorferl (|2005). 
All spectra were recorded with t he UV version of the Zurich 
Imaging Polarimeter ZIMPOL II dGandorfer et al.ll2004 at the 
1 .5 m McMath-Pierce facility of the National Solar Observatory 
on Kitt Peak (Arizona, USA). The spectrograph slit was set par- 
allel to the solar limb at the distance corresponding to the cosine 
of the heliocentric angle fi — 0.1 + 0.02. Due to the telescope 
guiding and seeing, the value of fi was not completely stable, re- 
sulting in a small uncertainty of the exact limb distance. Because 
the scattering polarization is very sensitive to small variations in 
very close to the solar limb, /i behaves like a semi-free pa- 
rameter for our analysis within strict limits defined by the error 
bars. 

One of the main problems for our fitting procedure is the un- 
certainty of the zero point of the polarization scale in the obser- 
vations. To solve this problem, we selected several wavelengths 
as continuum reference points, which are far away from very 
strong lines and least affecte d by blends. Based on the theoreti- 
cal continuum polarization (Stenflo 2005), we shifted the polar- 
ization scale of the observations appropriately. The same wave- 
length points were used to renormalize the I/I c data. 

We have chosen three wavelength regions for our calcu- 
lations. These regions were not observed simultaneously and 



correspond, in principle, to different areas of the solar surface. 
Therefore measured parameters can vary from region to region. 

3.3. Region 1: diagnostics near 3866 A 

The first wavelength region 3865.9-3866.8 A contains four CN 
violet system doublets with N - 6,1 and 52 from the (1,1) vi- 
brational band and with N -69 from the (0,0) vibrational band. 
Because of the relatively large Lande factors (g^ ~ I IN), the 
lines with N - 6 and 7 are very sensitive to the magnetic field 
and reach the saturated Hanle effect already at 5-10 G, while the 
lines with N-52 and 69 are much less affected by the magnetic 
field (see Fig. [2j. Therefore, this region is ideal for applying the 
differential Hanle effect technique, in particular since it contains 
lines with small and large N quantum numbers from the same vi- 
brational band, as the relative strength of the lines from different 
bands can be influenced by other factors. Furthermore, all lines 
can be measured simultaneously in a single observation because 
they are spread over less than 1 A. 

Figure[3]shows the accuracy of the fitted spectrum illustrated 
by ^-contours within the parameter space. All four free param- 
eters are well constrained by the data, and their best values are 
listed in Table Q] In particular we find a turbulent magnetic field 
strength of 27+5 G. The magnetic field B and the scattering 
probability c cs exhibit a somewhat linear dependence (Fig. [3] 
leftmost panel), which is easily understood. A greater scatter- 
ing probability leads to additional line polarization that has to be 
compensated with a stronger Hanle effect and thus a larger mag- 
netic field. Nonetheless, both B and c cs remain restricted to well 
defined limits because of the differential sensitivity to the Hanle 
effect of the involved CN lines and the roughly line independent 
scattering probability. Similarly, the parameters 7~~ nol - m and T e g 
also show a tendency towards a linear dependence, because a 
greater effective temperature reduces the line depth in the inten- 
sity spectrum that can be balanced by additional line opacity. 

The modeled spectra are compared with observations in 
Fig.|4] For the calculated curves we have assumed the best values 
of the parameters r e jf, c cs , and 7~ n0lm as listed in Table Q] while 
varying the magnetic field strength. The differential Hanle effect 
is clearly visible and constrains the magnetic field strength to 
well defined limits. It has to be stronger than 15 G to saturate 
the highly sensitive N -1 doublet. The upper limit for the field 
strength is defined by the only modest observed depolarizations 
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Fig. 5. x 2 - contours as in Fig. |4]but for the second wavelength region close to the (1,1) bandhead. We obtain a minimum value 
\L = 3.62. 
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Fig. 6. Fits to the observations (thin solid lines) of Stokes I/I c and Q/I c in the second wavelength region close to the (1,1) bandhead. 
The assumed magnetic field strengths in the theoretical spectra are G (no depolarization, dashed line), 25 G (best fit, thick solid 
line) and 45 G (too much depolarization, dotted line). Labels at the bottom of the lower panel indicate the N quantum numbers of 
the lower states and the line positions within each doublet (close to the bandhead N is omitted due to overlap). The labels at the top 
of the lower panel indicate the wavelength points which were selected for the fitting procedure. 



of the lines with N = 52 and 69. The observed doubling of the 
N = 52 doublet probably appears due to a strong Ti I blend ne- 
glected in our model. We obtain the best fit for the magnetic field 
strength of 27 G. 



3.4. Region 2: diagnostics in the (1,1) bandhead 

Molecular bandheads have a high potential for successful differ- 
ential Hanle effect diagnostics. They contain a mixture of lines 
with different /-numbers and correspondingly different mag- 
netic field sensitivities. On the other hand, all these lines belong 



to the same rotational branch and vibrational band so that many 
of the relevant line formation parameters nearly coincide. Here 
we present the results of our fitting procedure in the (1,1) band- 
head of the CN violet system. 

The fitting procedure is illustrated in Fig. [5] by the x 2 ' 
contours. Due to numerous neglected blends from atoms and 
molecules other than CN the quality of the fit is worse than in the 
first wavelength region (Sect. |3.3l l. Furthermore, the noise level 
in this region is also a bit higher. Nonetheless, we determine a 
magnetic field strength of 25+5 G in good agreement with the 
first region. 
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Fig. 8. Fits to the observations (thin solid lines) of the CN violet system in the third wavelength region close to the (0,0) bandhead. 
The magnetic field strengths in the modeled spectra were set to G (no depolarization, dashed line), 1 1 G (best fit, thick solid 
line), and 25 G (too much depolarization, dotted line). Labels at the bottom of the lower panel indicate the CN doublets with the N 
quantum numbers of the lower states. The labels at the top of the lower panel indicate the wavelength points chosen for the fitting 
procedure. 



The modeled spectra are compared to observations in Fig. [6] 
The bandhead of the (1,1) band is clearly visible near 3871.5 A. 
The selected wavelength region also contains a few lines from 
the (0,0) band. Contrary to the first region, lines with different 
Hanle effect sensitivities can not be easily distinguished as all 
lines are blended with each other. Even when neglecting the 
magnetic field it is possible to fit the overall average degree 
of polarization. However, the exact shape of the Q/I c spectrum 
can only be reproduced when accounting for a magnetic field 
strength within well defined limits. 

In this wavelength region we had to adopt a cosine of the 
heliocentric angle /i = 0.12, which still lies within the error limits 



of the observed data. In fact, we failed to obtain good fits in this 
region with fi = 0.10, because the continuum polarization in the 
incident radiation field was too large. 



3.5. Region 3: diagnostics close to the (0,0) bandhead 

The CN lines near the (0,0) bandhead are very strong and heav- 
ily blended. This results in a significant optical thickness larger 
than 1 at the (0,0) bandhead so that the single scattering approx- 
imation loses its validity. Therefore, we selected a region which 
is 4 A away from the bandhead. 



A. I. Shapiro et al.: Hanle effect in the CN violet system with LTE modeling 



9 



Table 1. Least squares fits to the magnetic field strength B, 
the optical depth normalization coefficient 7" norm , the scattering 
probability c sc , the temperature of the isothermal CN layer T e ff 
and cosine of the heliocentic angle yU in three considered regions 
of spectrum 



Region 


B,[G] 




c cs 


T c ft, K 




1 


27 ±5 


2.0 


0.375 


4620 


0.1 


2 


25 + 5 


3.5 


0.240 


4620 


0.12 


3 


11 ±2 


2.3 


0.142 


4460 


0.11 



This third wavelength region contains a mixture of lines with 
N ~ 10 and N ~ 40 from the P branch. The main features rele- 
vant for diagnosing the magnetic field in this spectral window 
are groups of neighboring lines containing a doublet with 10 
and one with N ~ 40, which are particularly well seen in the in- 
tensity spectrum. The doublets within one group can however be 
distinguished due to a consistent wavelength shift and different 
polarization characteristics. The lines with N ~10 have a larger 
effective polarizability W2 ~ 0.5 (compared to W2 close to the 
asymptotic limit of 0.1 for lines with A^~40), which in principle 
causes a larger line polarization in a non-magnetic medium. The 
Hanle effect introduces a competing contribution to the relative 
line polarization because the lines with small N are much more 
sensitive to the magnetic field. From Fig. [8] we see that the con- 
sistent modeling leads to smaller polarization and the expected 
stronger Hanle effect in the lines with small N. Therefore, we 
can indeed infer the magnetic field strength from the Q/I c spec- 
trum. 

The results of the fitting procedure are presented in Figs. [7] 
and [8] Again all the free parameters are well constrained and 
we reach a reasonable quality of the fit, except in three narrow 
regions due to missing blend lines in our model at 3879.0 A, 
3879.85 A, and 3880.2 A, which is clearly seen in the intensity 
spectrum. The best fit is obtained for a magnetic field strength 
of 11+2 G, which is significantly smaller than in the first two 
analyzed spectral windows. 

3.6. Discussion 

The above comparison of the modeled spectra to the observed 
data was performed independently in the three selected wave- 
length regions. The obtained sets of free parameters for each re- 
gion are collected in Table Q] They agree quite well in spite of 
the simple modeling, but there exist some differences that should 
be addressed. 

The magnetic field strength found in the third region lies sig- 
nificantly below the value determined for the first two spectral 
regions. This is however not necessarily a contradiction, because 
the analyzed data were not observed simultaneously and sam- 
pled in principle different spatial positions on the Sun. 

Furthermore, the line strength and accordingly the optical 
thickness changes from the first region to the third one. This in- 
fluences the formation heights of the CN lines, which leads to 
different physical conditions in the corresponding atmospheric 
layers. Lines in the third region are stronger than in the first two 
and therefore form higher in the atmosphere where the tempera- 
ture is smaller. A comparison of the contribution functions cal- 
culated for the typical lines in each spectral window shows that 
the difference in temperatures is around 100 K, consistent with 
the difference of 160 K obtained with our least-squares method 
(see Table[TJi. Our results suggest that the magnetic field strength 
may decrease with height, but a detailed interpretation of these 



small differences lies outside the scope of the current model and 
will be clarified by future work. 

In addition, the continuum normalization and polarization in 
the third region remains an issue due to broad wings of nearby 
strong atomic lines including the H8 line from the Balmer se- 
ries. The intensity normalization has little effect on the inferred 
magnetic field strength because it affects both the theoretical and 
observed Q/I c signal in the same way. However, the additional 
opacity from those nearby atomic lines would probably reduce 
the continuum polarization and slightly increase the magnetic 
field strength by a few Gauss. 

Overall, the physical parameters determined with the fit lie 
in a reasonable range. The small scattering probability is con- 
sistent with expectations and with a relatively small influence of 
non-LTE effects. The total optical thickness of the CN layer is 
close to but still smaller than 1 (maximal optical thicknesses in 
the considered regions are 0.26, 0.5 and 0.32 for the first, second 
and third regions correspondingly) in accordance with our as- 
sumptions. Finally, the temperature T e s of the CN layer around 
4500-4600 K that is required for the best fits, corresponds well 
to the temperature in the upper solar photosphere where the CN 
lines are formed. 

An interesting extension of our method is the interpre- 
tation of the center-to-limb variations of the continuum and 
line radiation field, as it was done for Sr I 4607 A line by 
IStenflo et aT] d 19971) and for t he C 2 and MgH molecular lines 
bv lFaurobert & Arnaudl(l2003l) . Such observation for the CN vi- 
olet system are not available yet, but could be very useful for 
constraining the atmospheric models. 



4. Conclusions 

We have employed the differential Hanle effect for the interpre- 
tation of observations of the CN violet system, which allowed us 
to determine the turbulent magnetic field strength in the upper 
solar photosphere. This expands the number of molecules used 
for Hanle effect diagnostics, which was so far only based on C2 
and MgH. The main conclusions can be summarized as follows: 

- We have developed a new model that can diagnose the tur- 
bulent magnetic field via the differential Hanle effect. It is 
applicable for optically thin lines, relies on simple and few 
assumptions, and consistently accounts for continuum polar- 
ization, depolarization of the continuum by line opacity, in- 
trinsic line polarization by coherent scattering, and the Hanle 
effect. 

- Lines from the CN violet system are well suited as a diagnos- 
tic tool based on the differential Hanle effect. They are very 
prominent in the second solar spectrum, and there exist sev- 
eral narrow spectral regions containing simultaneously lines 
with small and large J quantum numbers. These lines exhibit 
different sensitivities to the Hanle effect, since the effective 
Lande factor is inversely proportional to J in the CN violet 
system. 

- We have identified three spectral windows that are ideal 
for Hanle effect diagnostics and contain a small number of 
atomic and molecular blends. 

- In all three spectral regions we obtained a good quality of 
the fit and determined turbulent magnetic field strengths of 
27+5 G, 25±5 G, and 1 1+2 G, respectively. 
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Table 2. Line parameters used for the fitting procedure in the first spectral region (3865.9-3866.8 A). The columns list the wave- 
length (A), the line identification (ID), the polarizability (W2), the lifetime of the upper level (tr), the turbulent magnetic field 
strength at which the line polarization is reduced to 30% of its non-magnetic value neglecting elastic collision (B30), the vibrational 
number which coincides for the upper and lower states (v), and the oscillator strength (f v js). 



A, A 


ID 


w 2 


t r , 10~ 8 s 


P30 


V 


JvJS 


A, A 


ID 


w 2 


Tr, 10~ 8 S 


P30 


V 


fvJS 


3865.978 


«i(12) 


0.1239 


7.2719 


22 





0.391 


3865.998 


#2(12) 


0.1260 


7.2736 


22 





0.362 


3866.013 


2i2(12) 


0.3982 


7.2736 


22 





0.001 


3866.096 


Pi (6) 


0.0604 


7.4397 


9 


1 


0.154 


3866.090 


621(6) 


0.3916 


7.4397 


9 


1 


0.002 


3866.098 


Pi(6) 


0.0545 


7.4635 


9 


1 


0.130 


3866.160 


Pi (69) 


0.0957 


7.0063 


111 





2.016 


3866.086 


221 (69) 


0.3999 


7.0063 


111 





0.000 


3866.235 


P 2 (69) 


0.0957 


7.0069 


111 





1.987 


3866.400 


Pi (52) 


0.0944 


7.2450 


86 


1 


1.254 


3866.344 


221 (52) 


0.3999 


7.2450 


86 


1 


0.000 


3866.460 


P 2 (52) 


0.0943 


7.2456 


86 


1 


1.230 


3866.540 


Pi(7) 


0.0650 


7.4251 


11 


1 


0.178 


3866.533 


22! (7) 


0.3938 


7.4251 


11 


1 


0.002 



Table 3. Line parameters used for the fitting procedure in the second spectral region (in the (1,1) bandhead), with the same param- 
eters as in Table [2] 



A, A 


ID 


W 2 


tr, 10~ 8 s 


P30 


V 


fvJS 


A, A 


ID 


W 2 


Tr, 10- 8 S 


P30 


V fvJS 


3869.672 


Pi(16) 


0.0829 


7.3832 


26 


1 


0.394 


3869.655 


22,(16) 


0.3988 


7.3832 


26 


1 0.001 


3869.695 


P 2 (16) 


0.0819 


7.3845 


26 


1 


0.370 


3869.831 


Pi (42) 


0.0931 


7.2954 


70 


1 1.015 


3869.786 


221 (42) 


0.3998 


7.2954 


70 


1 


0.000 


3869.886 


P,(42) 


0.0929 


7.2959 


70 


1 0.991 


3869.920 


P-i(7) 


0.1397 


7.2917 


14 





0.246 


3869.938 


R 2 (T> 


0.1457 


7.2980 


14 


0.217 


3869.947 


2l2(7) 


0.3953 


7.2980 


14 





0.002 


3869.916 


P,(17) 


0.0838 


7.3806 


28 


1 0.418 


3869.898 


22l(17) 


0.3990 


7.3806 


28 


1 


0.001 


3869.940 


P 2 (17) 


0.0829 


7.3817 


28 


1 0.394 


3870.070 


Pi (64) 


0.0954 


7.0428 


103 





1.871 


3870.001 


22,(64) 


0.3999 


7.0428 


103 


0.000 


3870.146 


P 2 (64) 


0.0953 


7.0433 


103 





1.842 


3870.070 


Pi(41) 


0.0929 


7.2999 


68 


1 0.991 


3870.026 


22,(41) 


0.3998 


7.2999 


68 


1 


0.000 


3870.118 


P 2 (41) 


0.0928 


7.3004 


68 


1 0.967 


3870.146 


P,(18) 


0.0846 


7.3780 


29 


1 


0.442 


3870.127 


22l(18) 


0.3991 


7.3780 


29 


1 0.001 


3870.176 


P 2 (18) 


0.0838 


7.3790 


29 


1 


0.418 


3870.282 


Pi (40) 


0.0928 


7.3043 


66 


1 0.967 


3870.239 


221 (40) 


0.3998 


7.3043 


66 


1 


0.000 


3870.340 


/M40) 


0.0926 


7.3048 


66 


1 0.944 


3870.460 


Pi(19) 


0.0854 


7.3754 


31 


1 


0.466 


3870.440 


22ld9) 


0.3992 


7.3754 


31 


1 0.001 


3870.489 


P 2 (19) 


0.0846 


7.3763 


31 


1 


0.442 


3870.530 


Pi (39) 


0.0926 


7.3086 


65 


1 0.944 


3870.488 


221 (39) 


0.3998 


7.3086 


65 


1 


0.000 


3870.580 


P 2 (39) 


0.0924 


7.3091 


65 


1 0.920 


3870.551 


Pi (20) 


0.0861 


7.3728 


33 


1 


0.489 


3870.530 


221 (20) 


0.3992 


7.3728 


33 


1 0.001 


3870.582 


P 2 (20) 


0.0854 


7.3736 


33 


1 


0.466 


3870.657 


Pi(6) 


0.1457 


7.2987 


12 


0.217 


3870.664 


P 2 (6) 


0.1538 


7.3078 


12 





0.187 


3870.672 


2n(6) 


0.3938 


7.3078 


12 


0.002 


3870.715 


Pi(21) 


0.0867 


7.3701 


35 


1 


0.513 


3870.693 


2 2 i(21) 


0.3993 


7.3701 


35 


1 0.001 


3870.751 


P 2 (21) 


0.0861 


7.3708 


35 


1 


0.489 


3870.662 


Pi (38) 


0.0924 


7.3128 


63 


1 0.920 


3870.621 


22,(38) 


0.3998 


7.3128 


63 


1 


0.000 


3870.715 


P2O8) 


0.0922 


7.3133 


63 


1 0.896 


3870.797 


Pi (63) 


0.0953 


7.0498 


102 





1.842 


3870.729 


22i(63) 


0.3999 


7.0498 


102 


0.000 


3870.871 


P 2 (63) 


0.0953 


7.0503 


102 





1.813 


3870.872 


Pi (22) 


0.0872 


7.3674 


36 


1 0.537 


3870.849 


22,(22) 


0.3994 


7.3674 


36 


1 


0.001 


3870.906 


/M22) 


0.0867 


7.3681 


36 


1 0.513 


3870.817 


P,(37) 


0.0922 


7.3170 


61 


1 


0.896 


3870.777 


22i(37) 


0.3998 


7.3170 


61 


1 0.000 


3870.870 


P 2 (37) 


0.0920 


7.3174 


61 


1 


0.872 


3871.007 


Pi (23) 


0.0878 


7.3647 


38 


1 0.561 


3870.983 


22,(23) 


0.3994 


7.3647 


38 


1 


0.001 


3871.042 


P 2 (23) 


0.0872 


7.3653 


38 


1 0.537 


3870.956 


Pi (36) 


0.0920 


7.3210 


60 


1 


0.872 


3870.917 


221 (36) 


0.3998 


7.3210 


60 


1 0.000 


3871.107 


P2O6) 


0.0918 


7.3214 


60 


1 


0.848 


3871.125 


Pi (24) 


0.0882 


7.3618 


40 


1 0.585 


3871.099 


22,(24) 


0.3995 


7.3618 


40 


1 


0.000 


3871.159 


Po(24) 


0.0878 


7.3624 


40 


1 0.561 


3871.222 


Pi (25) 


0.0887 


7.3589 


41 


1 


0.609 


3871.195 


221 (25) 


0.3995 


7.3589 


41 


1 0.000 


3871.079 


P,(35) 


0.0918 


7.3249 


58 


1 


0.848 


3871.042 


22i(35) 


0.3998 


7.3249 


58 


1 0.000 


3871.258 


P 2 (25) 


0.0882 


7.3594 


41 


1 


0.585 


3871.298 


Pi (26) 


0.0891 


7.3559 


43 


1 0.633 


3871.270 


22,(26) 


0.3996 


7.3559 


43 


1 


0.000 


3871.125 


P 2 (35) 


0.0916 


7.3254 


58 


1 0.824 


3871.382 


P 2 (27) 


0.0891 


7.3533 


45 


1 


0.633 


3871.417 


Po(28) 


0.0895 


7.3501 


46 


1 0.657 


3871.338 


Pi (27) 


0.0895 


7.3528 


45 


1 


0.657 


3871.309 


221 (27) 


0.3996 


7.3528 


45 


1 0.000 


3871.182 


Pi (34) 


0.0916 


7.3287 


56 


1 


0.824 


3871.145 


221 (34) 


0.3997 


7.3287 


56 


1 0.000 


3871.329 


P 2 (26) 


0.0887 


7.3564 


43 


1 


0.609 


3871.364 


Pi(5) 


0.1538 


7.3083 


11 


0.187 


3871.372 


P 2 (5) 


0.1655 


7.3224 


11 





0.158 


3871.379 


2l2(5) 


0.3916 


7.3224 


11 


0.002 


3871.222 


P 2 (34) 


0.0913 


7.3292 


56 


1 


0.800 


3871.258 


Pi (33) 


0.0913 


7.3324 


55 


1 0.800 


3871.223 


22,(33) 


0.3997 


7.3324 


55 


1 


0.000 


3871.298 


P 2 (33) 


0.0910 


7.3329 


55 


1 0.776 


3871.393 


Pi (30) 


0.0905 


7.3430 


50 


1 


0.728 


3871.361 


22i(30) 


0.3997 


7.3430 


50 


1 0.000 


3871.438 


P 2 (30) 


0.0902 


7.3435 


50 


1 


0.705 


3871.393 


Pi (29) 


0.0902 


7.3464 


48 


1 0.705 


3871.362 


22,(29) 


0.3996 


7.3464 


48 


1 


0.000 


3871.438 


P 2 (29) 


0.0898 


7.3469 


48 


1 0.681 


3871.382 


Pi (28) 


0.0898 


7.3496 


46 


1 


0.681 


3871.352 


221 (28) 


0.3996 


7.3496 


46 


1 0.000 


3871.366 


Pi(31) 


0.0908 


7.3396 


51 


1 


0.752 


3871.333 


2 2 i(31) 


0.3997 


7.3396 


51 


1 0.000 


3871.417 


P 2 (31) 


0.0905 


7.3401 


51 


1 


0.728 


3871.329 


Pi (32) 


0.0910 


7.3361 


53 


1 0.776 


3871.295 


22,(32) 


0.3997 


7.3361 


53 


1 


0.000 


3871.360 


/M32) 


0.0908 


7.3365 


53 


1 0.752 
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Table 4. Line parameters used for the fitting procedure in the third spectral region (close to the (0,0) bandhead), with the same 
parameters as in Table [2] 



A, A 


ID 


Wi 


t r , 10~ 8 s 


P30 


V 


fvJS 


A, A 


ID 


W 2 


Tr, 10~ 8 S 


P30 


V 


fvJS 


3879.174 


p 2 m 


0.0650 


7.3078 


12 





0.217 


3879.185 


Pi (8) 


0.0686 


7.2987 


12 





0.246 


3879.177 


Q 2 t(8) 


0.3953 


7.2987 


12 





0.002 


3879.247 


Pi (48) 


0.0939 


7.1431 


78 





1.407 


3879.195 


22i (48) 


0.3999 


7.1431 


78 





0.000 


3879.306 


P 2 (48) 


0.0938 


7.1435 


78 





1.378 


3879.566 


Pi(9) 


0.0716 


7.2917 


14 





0.275 


3879.557 


221 (9) 


0.3963 


7.2917 


14 





0.002 


3879.580 


P 2 (9) 


0.0686 


7.2980 


14 





0.246 


3879.649 


Pi (47) 


0.0938 


7.1485 


77 





1.378 


3879.598 


g 2 i (47) 


0.3999 


7.1485 


77 





0.000 


3879.707 


P 2 (47) 


0.0937 


7.1489 


77 





1.349 


3879.952 


Pi(10) 


0.0740 


7.2863 


16 





0.304 


3879.942 


22l(10) 


0.3970 


7.2863 


16 





0.001 


3879.967 


Pi(lO) 


0.0716 


7.2909 


16 





0.275 


3880.029 


Pi (46) 


0.0937 


7.1538 


75 





1.349 


3879.980 


22i (46) 


0.3999 


7.1538 


75 





0.000 


3880.087 


P. (46) 


0.0935 


7.1542 


75 





1.320 


3880.317 


/Mil) 


0.0761 


7.2819 


17 





0.333 


3880.305 


221(H) 


0.3975 


7.2819 


17 





0.001 


3880.332 


P 2 (ll) 


0.0740 


7.2854 


17 





0.304 


3880.389 


Pi (45) 


0.0935 


7.1590 


73 





1.320 


3880.340 


22i (45) 


0.3999 


7.1590 


73 





0.000 


3880.446 


P 2 (45) 


0.0934 


7.1594 


73 





1.291 


3880.662 


Pn(12) 


0.0778 


7.2782 


19 





0.362 


3880.649 


22l(12) 


0.3979 


7.2782 


19 





0.001 


3880.678 


P 2 (12) 


0.0761 


7.2809 


19 





0.333 


3880.728 


Pi (44) 


0.0934 


7.1641 


72 





1.291 


3880.681 


221 (44) 


0.3998 


7.1641 


72 





0.000 


3880.784 


P,(44) 


0.0933 


7.1645 


72 





1.262 


3880.985 


Pi(13) 


0.0794 


7.2749 


21 





0.391 


3880.971 


22.(13) 


0.3982 


7.2749 


21 





0.001 



